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HAT I DUAL ADVISORY COMMITTEE IP OR AERONAUTICS 
TECHNICAL NOTE NO. 970 

ON THE CIRCULATORY SUBSONIC EL ON OE a COMPRESSIBLE ELUID 

PAST A CIRCULAR CYLINDER 
By Lipman Bers 

* 

SUMMARY 


The circulatory subsonic flow around an infinite circu- 
lar cylinder is computed using the linearized pres sure -volume 
relation, hy a method developed in a previous report. Eor- 
mulas and graphs are given for the velocity and pressure dis- 
tributions, the circulation, the lift, and the dependence of 
the critical Mach number upon the position of the stagnation 
p oint . 


INTRODUCTION 


The rigorous solution of the differential equations of 
a two-dimensional steady potential compressible gas flow in- 
volves considerable mathematical difficulties. Numerical in- 
tegration is rather laborious and hardly ever yields results 
of a general character. Therefore, considerable attention 
was given to approximate analytical methods. 

Tchaplygin (reference l) noticed that the differential 
equation of the velocity potential takes a rather simple form 
if the exponent Y in the polytropic pressure-density rela- 
tion is replaced by -1. The equation then becomes the well- 
known equation of a minimal surface. Furthermore, in the so- 
called hodograph plane, the equation may be transformed into 
the Laplace equation. The physical meaning of the linearized 
pressure-volume relation (Y = -l) has been clarified by 
Busemann (reference 2} and especially by Von Karman (refer- 
ence 3) and Tsien (reference 4). They showed that using this 
relation amounts to replacing the actual pressure-volume curve 
by its tangent. Another way of justifying the use of the lin- 
earized equation of state is indicated in this report. 
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Tsien derived a formula transforming a given incompress 
ible flow around a closed profile into a compressible flow 
(satisfying the linearized equation of state; around another 
(slightly distorted) profile. However, this formula is ap- 
plicable only to circulation-free flows. In a previous re- 
port (reference 5) the author developed formulas which per- 
mit the construction of circulatory flows as well. It has 
been shown that under very general conditions every incom- 
pressible flow yields a compressible flow and that all com- 
pressible flows may be obtained in this way. 

In the present report, this method is tested and illus- 
trated by computing the circulatory flow around a circular 
cylinder. (The cir culat i on-f r ee flow around a circular cyl- 
inder has been treated by Tamada, who used Tsien’ s method.) 
(See reference 6.5 The circular profile, although of no 
technical importance, has been chosen because of the simplic 
ity of the computations and the possibility of comparing re- 
sults with those obtained by different methods. 

It should be emphasized that the use of the linearized 
pressure-volume relation restricts the investigation to 
purely subsonic flows. 

The investigation, conducted at Brown University, was 
sponsored by and carried out with the financial assistance 
of the National Advisory Committee for Aeronautics. 

The author is indebted to Mr. Charles Saltzer for val- 
uable assistance. 


SYMBOLS 


a local speed of sound 

a Q speed of sound at a stagnation point 
C L lift coefficient 

C( £) normalized complex potential of an incompressible flow 
in the £-plane 

L lift per unit span 

M 


local Mach number 
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M stream Mach number (Mach number at infinity) 

CD 

M^c critical stream Mach number 

n parameter defined by the equation n = J\ +. (qja 0 ) s 

P profile in thd z-plane (the plane of the compressible 
flow) 

p pressure 

p Q pressure at the stagnation point 

q local speed of a compressible flow 
q g value of q for which q = a 

a speed of the compressible flow at Infinity 

£ radius of the circular cylinder 

Hi auxiliary parameter defined by the equation 
E x = cos S/cos P 

Ep radius of curvature of the profile P 

Ep- radius of curvature of the profile T 

r parameter defined by the equation 2r/(l - r S ) = q s /a 0 
S pressure coefficient 

u,v components of the local velocity of the compressible flow 

z = x + iy complex variable in the plane of the compressible 

flow 

Z lt Z 2 ,Z 3 auxiliary complex -.ariables 

a absolute value of the argument of a stagnation point of 
the compressible flow about a circular cylinder 
(angle of attack) 

P,5iP 1 ,P 2 angles used to locate stagnation points of auxil- 
iary flows In the £ -plane 
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y exponent in the polytropic relation 

T c circulation of the compressible flow 

Ti circulation of the incompressible flow 

g angle defined by the equation 8 = tt/2 - £ 

£ = £ + iT| complex variable in the plane of the auxiliary 
incompressible flow 

t 

A argument of a point on the circle iZJ = R 
(A parameter defined by the equation n 3 = (n - l)/(n + l) 

fl profile in the plane of the incompressible flow (£-plane) 

p density 

p Q stagnation density 

cr absolute value of the argument of a stagnation point of 
the flow in the Z-plane 

vj/ angles used to locate points on the profile n 

cp velocity potential of the compressible flow 

"X, angle defined by the equation X= (tt/ 2 - fi)/('2 - n) 

oj argument of a point on the cross section of the circular 

cylinder 


ANALYSIS 

1. Reduction of the Problem to the Differential Equation 

of a Minimal Surface 


Let u and v denote the x and y components of the 
velocity of a steady two-dimensional flow of a compressible 
fluid. If the flow is irr ot at i onal , a velocity potential 
cp(x f y) may be introduced, such that 


o 


$ 3 2 
3x 


a. 


dco 


u = a 


v 


( 1 ) 
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/ 

where a 0 is the speed of sound at a stagnation point. The 
continuity equation yields the equation 


0 


3x V r 3x/ 3y V r 3y/ 


0 


0 


(3) 


p being the density. 

If the flow occupies the domain exterior to a closed 
profile P and at infinit 3 r approaches a uniform flow in the 
positive x direction, <p must satisfy the boundary condi- 
tions 


3x a Q ’ ay 


0 as x® + y s 


eo 


(3) 


where q^ is the speed at infinity (speed of the undisturbed 
f low) , and 


I s = 0 on P 

031 


(4) 


(JL. denotes differentiation in the direction normal to the 
\3n v 

profile Pj, 

The speed of the fluid q is given by 


3 3 s 3 

q = u + v = a 0 


(MT + 


and the Kutta- Joukowski condition requires that 

q < oo 


(5) 


It will be assumed that the pressure p and the density 
are connected by the relation 


P 


A + Bp 


Y 


( 6 ) 


where A, B, and Y are constants. Por all physically im- 
portant cases 


1 < Y < 2 
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The standard value of 7 for air is 1.405. Equation (6) and 
Bernoulli's theorem imply that 



(7) 


S 8 V *“ 1 _ 2 

a “ a o « 1 


M 3 = 


(q./ a 0 )' 


i - 


7-1 


U/ a o) 


( 8 ) 

(9) 


Here p 0 is the density at a stagnation point, a = dp/dp, 

the local speed of sound and M = <i/a, the local Mach number 
Tt is seen that the flow is subsonic as long as 


q> < q, 


= a 


/- ' .. a — 

V 7 + 1 


( 10 ) 


Equation (7) can be written as a relation (depending on 
the parameter 7) between the two dimensionless quantities 

p/p c and <l/ a 0 5 



(11) 


Equation (s) can now be written as 



( 2 1 ) 


Expanding the right side of (ll) into a power 
q/a Q yields 


series in 
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f(x, ■»).!. ip-Y ♦ j 

^a 0 / 2\a 0 / 8. \a 0 / 


— + . 


Evidently, for small values of q./a Q the relation between p 
and q./a Q is only slightly affected by the value of V. In 
particular, if the actual value of V is replaced by -1 , the 
resulting error in p/p 0 is of the same order of magnitude 
as 


y + i 
8 



Table I gives the values of p/p 0 for ^ = 1.405 and 
Y = -1 and for values of q/a 0 from 0 to 0.912 (0.912 is 
the value of q s /a Q for Y = 1.405). 


On the basis of the foregoing remarks it may be con- 
cluded that for flows where q/a Q is Bmall only a slight 
error will be made if the actual value of Y is replaced by 
-1 in (ll). Since 


(fe’ 



1 

yr + ( q./a 0 ) 3 


TABLE I.- DENSITY AS A FUNCTION OF SPEED 


Q . 

a o 

0 


0.2 

0.3 

0.4 



0.5 


0.7 

O.S 

0.9 

0.912 

y = 1.405 

q. 

PS 

1.000 

0.995 

0.920 

0.956 

0.922 

0.380 

0.330 



0-773 

0.710 

0.643 

0.634 

y = -i.ooo 

E 

1.000 

0.995 

0.980 

0.95s 

0.929! 0 .S 94 

i 

0 .S 5 S 

0 .S 20 

0.7S0 

0.743 

0.739 


equation (2 1 ) becomes 
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(l + C^-f\ §fse. - 2 §S£. §2. §LSE- +(l + = 0 (12) 

v. i 3i 2 Bx By BxBy (• ''Bx' J By 3 

This is the well-}cnown equation of a minimal surface, and its 
solution can he expressed in terms of analytic functions of 
a complex variable. 

After the values of q/a 0 have been found by integrat- 
ing the approximate equation (12) under the boundary condi- 
tions (3), (4), and ( 5 ), the values of the density and the 
looal Mach number must be determined by means of formulas 
(7) and (9) with the correct value of V. 

The use of the approximate equation (12) can also be 
justified by the remark that the constants A and B can be 
determined so that the curve (in the p, p-plane) given by 
the linear pressure-volume relation 

p = A + — (13) 

P 

will be at some -point tangent to the curve given by the actual 
pressure-density relation 



(see references 3 and 4). The linear pressure-volume relation 
was introduced by Tchaplygin and has been used by Eusemann, 
Demtehenko, Von Karman, and Tsien. 


2. Formulas for the Solution of the Approximate Problem 

In a previous report (see reference 5) a method was 
given which permits the construction of solutions of equation 
(12) satisfying the boundary conditions (3), (4), and (5). 

It has been shown that this method yields solutions for all 
profiles P possessing at most two sharp edges. The only 
restrictive condition concerns the value of q.eo/ a o which 

should not exceed ,*/■— - 0.866. Now for V * 1.405, = 0,912. 

Thus, the method would fail only for profiles for which the 
maximum local speed is very close to the speed of the undis- 
turbed flow - that is, for very thin profiles. 
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This method is based upon a transformation of an incom- 
pressible flow into a solution of equation (12). 

let n be a closed profile in the £-plane, (£ = | + iTt) 
possessing at most two edges. Let &(£ ) be a normalized com- 
plex potential (i.e., the function G- = [velocity potential 
+ i( stream function) j such that the speed at infinity equals 
l) of an incompressible flow around the profile H ; G(£ ) is 
an analytic function defined in the domain exterior to rt and 
satisfies the conditions 

&*(“) = 1 

Im G- = constant on n 


as well as the Kutta- Joukowsky condition 


. 10*1 


max 


< ® 


(Im = imaginary part of). Let n be a real number such that 

1 < n < 2 

The function 

(2.1) z = x + iy = — 


±Ji o*a) x ” tt d£ - a£ (14) 


(the bar denotes the conjugate complex quantity) maps the do- 
main exterior to n into a domain exterior to a closed pro- 
file P. The function 


cp = «/n 2 - 


1 Re G 


(15) 


(Re = real part of) considered. as a function of x and y 
satisfies the differential equation (12) and the boundary 
conditions, (4), (5), and (6), where 


Is 

a 0 


- 1 


(16) 


Thus cp may be considered as the potential of a compressible 
flow around P. The speed q of the compressible flow (at a 
point z corresponding to a point £) is given by 
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q. 

a o 


2 (rr t)* i°'t» ^ 

x - fs-^_i v )1 <; ’ { £)l s/n 

Vn + 1/ 

3 ^/a 0 i&«(Ql. l/n 
n + 1 - (n - 1) JO' (O 2 ^ 11 


(17) 


The angle "between the direction of the compressible flow and 
the i-axis, 

8 = tan -1 — 
u 

is given by 


6 = — arg &» ( £) 


(18) 


(arg = argument 
P is l/n t ime s 
ing point . If 
a point z and 
sponding point 


of). In particular, the slope of the profile 
the slope of the profile n , at a correspond- 
Rp denotes the radius of curvature of P at 
Rp) the radius of curvature at the corre- 
£ of n , then 




The profile P can be constructed graphically using this 
informat i on . 

The proofs of these assertions will be found in the re- 
port quoted above. 

Hote that for values of n which are close to 1 (i.e., 
for small values of q cc / a o^ ik® profile P will be slightly 

different from n . Por, then, the coefficient of the first 
integral in (14) is close to 1, the exponent of the integrand 
is close to 0, so that the first term on the right side of 
(14) is close to £, and the second term is small compared to 
the first. 


However, the slopes of P and n (at corresponding 
points) are different. If it is desired that the profile P 
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should possess at the two stagnation points angles £ 3 . and 
0 2 , the profile n must possess at its stagnation points 
* the angles n 0 x and n 0 2 . 

It will he noted that the method does not yet permit the 
solution of the boundary value problem for a given profile. 
However, it is possible to choose the "conjugate profile" n 
in such a manner that the profile P will differ very little 
from a given profile. 

Remark ; If G- is a one-valued function (i.e., for a 
flow without circulation), the formulas (14), (15), (17), 

(18), and (19) can be replaced by the simpler relations which 
are equivalent to the formulas found by Tsien (see reference 


4) 

z = ~~ t - f ( 30) 

cp = J n S - 1 He G (21) 

2 ^T IS ' ( ^' (22) 

1 - (rtr T 1 s,( S )|3 

8 = tan" 1 &' (£) (23) 

a p - 1 ( t,) t a j R n (24) 


However, if there is circulation, formula (20) does not yield 
a closed profile. 

3. Incompressible Plow Yielding a Compressible Plow 
around a Nearly Circular Cylinder 


The method described In the foregoing will be applied 
presently to the construction of a circulatory compressible 
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flow around a nearly circular profile P. It is natural to 
start with a profile n consisting of an upper arc of the 
unit circle j£l = 1 and a lower arc of a circle which in- 
tersects the first arc at the points £ = e**^® and £ = -e^® 
and forms the angles nir with the first arc. (See fig. 1.) 
The center of the lower arc Is situated at £ = i(R x sin £ 

- sin. 8.) where 


g = S-(n-l)tr 


(25) 


Its radius is equal to 


Rj. = cos 6/eos g 


It may "be assumed without loss of generality that 0 < 8 < ?■ 

2 

and Rj < 1 . 


Let Gr( £ ) he the complex potential of an incompre s sihle 
flow past n which possesses stagnation points at £ = 
and £ = and has the velocity 1 at infinity. 

In order to compute &(£) the domain exterior to n is 
mapped conformally into the domain \Z[ > R in an auxiliary 
Z-plane, hy a transformation which takes £ = oo into Z = ce 
and satisfies the condition dZ/d£ >0 at infinity. The 

points £ = e"*^ and £ = -e*^ are taken into the points 
Z - Re-* 0 " and Z = -Re^, respectively, where 

<r . — i (a?) 

2 - n 2 - n 2 

Set 

. R - (as) 

(2 - n)cos <7 

then (dZ/d£)^ =oj = 1. Hence &(£>, considered ae a function 
of Z, is given hy 

P 2 

+ — + 2 Ri sin or log Z 
Z 


G- = Z 


(29) 
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At a point £ = ® o f tie upper arc of n the speed 

of the incompressible fluid is given by 


G‘(e i<5 ) = 2B S < sia , A + , ein g l S 
j s i n $ + sin 6 1 


(30) 


At a point £ = iC^ sin 3 - sin 6) + R e* - ^ of the lover ar< 
of n the speed of the incompressible flow is given by 

lo’fiCRj^ sin 3 - sin S) + = 2 ~ ( 3] 

It R* sin y\r + sin p| 

In these formulas A is the argument of the point on the 
circle |Z| = R into which the point £ is taken. The an- 
gle -(A- cr) can be determined by the formulas 


tan —(a - cr) = 
2 


cos cr 


< 

[ 1 

COB - 

l 2 


($ - 6 


(for a point on the upper arc) and 


tan i( A - cr) = - 1 

2 co 8 a 


sin i(\{r + p ) ' 3-11 

cos i(ij/ - 3 ) 

2 


+ sin a W33) 


(for a point ofi the lower arc). The details of the computa- 
tion will be found in the appendix. 

The maximum of | G- 1 j is reached at the top of the pro- 
file n, where $ = A = II . Hence 

3 


]G-'i_ = 2(l + sin 8 )R 


cos 4 — (tt - 2<r) 
s 4 

COS 4— ( TT - 28) 
4 


This is easily transformed inti 


(2 - n) ‘ 


(l - sin 6)co 


t 2 nln ~ — 


2(2 - n) 


(34) 
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The incompressible flow around n is now transformed 
into a compressible flow around the profile P (in the z- 
plane). The transformation is given by the formulas of the 
preceding section. By virtue ef (l?) and the actual equation 
of state, the maximal local Mach number will be equal to 1 
if 

=C/f-H r )° 

where r is determined from the equation 


2 r , I 2 

* = q / a s / 

1 - r 3 s' o 


V + 1 


For V = 1.405 


r = 0.3875 . . . 

Thus, for a given S, the maximal admissible n will be 
given by the relation 



sln 6 )eot 3 = ( 2 - r n 

2(2 - n) V* - 1 / 


Set 


6 = | - 2 e ' (35) 

then e is to be determined as the positive root of the 
equat ion 

A(n) sin e = tan — (36) 

2 - n 

with 

A(n) = gT-tt/a ( n.-. l W 4 ( 37 ) 

2 - n ^n + l' 


Numerical values of n and e determined by this eauation 
(for V = 1.405) are given in table IX. It is seen that the 
values of n are very close to 1 . 
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The radius of curvature of the profile P can. be com- 
puted by means of formula (19). The result of this computa- 
tion. for varioiis values of n and 6 are given in table 

III. It will be noted that the radius of curvature of P 

is nearly constant. Therefore, the profile P is nearly a 

circle of unit, radius; A point £ = e*® of the upper arc 
Is taken appr oximat ely into the point e^ with 

(!>=—■$- - 1 - TL (38) 

n a 2 

A point £ = i(R^ s ^ n P - s ^ n 8) + of the lower circle 

is taken approximately into the point e* ,Ji with 


U) 


\£_ + n ~ 1 tt 
n n 2 


(39) 


This follows from the fact that the slope of P at a point 
z must equal l/n times the slupe of n at a corresponding 
point £. In particular, the stagnation points of the com- 
pressible flow will be situated at e”* 0, and -e ict ' l where 


6_ n - 1 tt 

n n 2 


(40) 


a will be positive for 6 > (n - 1)tt/2. 

Greater accuracy could be achieved by taking as the 
radius of P the arithmetic mean of the values of Rp , or 

to obtain this radius by a graphical construction of P. 
However, this correction seems to be too insignificant to 
justify the additional computational labor. Of course, the 
size of P is of no importance as far as the velocity dis- 
tribution is concerned. 

The graphical construction of the profile P Is shown 
in figur'e 2. Figure 3 contains several examples of profiles 
P constructed by this method. The deviation of P from a 
circle is remarkably small. 
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4. Critical Mach Number 

The critical stream Mach number M. is defined as 

the value of H a for which the maximal local Mach number 

eauals 1.. In order to compute the critical Mach number for 
a flow around a circle with the angle of attack a (the 
angle of attack being defined as the negative argument of a 
stagnation point), set 


a 


g - n tt 2e 

n 2 n 


(41) 


and 

M 2 * (42) 

1 - 1(V - 1 ) ( n 3 - 1) 

2 

where n and e are connected by equation (36). [Cf. (9), 
(16), and' (40).j In this way the values of , c as a 
function of a (for V = 1.406) given in table IV have been 
computed. The relation between ^a>,c ail< 2- ct is plotted in 
figure 4. 


5. Velocity and Pressure Distribution 

In order to compute the velocity distribution of a com- 
pressible flow along a circular profile, set 


l* * 1 


Jk. 


i + ^m; 


(43) 


6 = net, + (n - 1 )— 

2 


(44) 


where is the desired stream Mach number and a the de- 

sired angle of attack. Then the dimensionless speed of the 

compressible flow at a point e ^ of the circle is given by 
formula (17). And J G- * | must be computed by means of for- 
mulas (30) and (3l) . The first is to be used for 2 > u> > -a, 

Ct — * 
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the second for 


H < uq < -a. 
2 — ~ 


(Note that the velocity dis- 


tribution is symmetrical with respect to the vertical axis.) 
The values of $ and \}r corresponding to a given value of 
A are given by equations (38) and (39). In this way the 
velocity distribution plotted in figure 5 has been computed. 


Once the velocity distribution Is known, the pressure 
distribution may be determined by means of equations (7) and 
the actual equation of state. The dimensionless pressure 
coefficient 


is given by 


S = 


” P 


— p 
2 ° 


CO ^-03 


S = £- 


-(*- 

- 


1 - _i _aL2\ 

2 


Y/(Y~i) 


a 0 3/ 


v - i iqfYr^r 


an®/ 


(45) 


For an incompressible flow 


S 



(46) 


The values of S plotted in figure 6 have been computed by 
these formulas. 


6. Circulation and Lift 


flow 


By 

i s 


virtue of (l) the circulation 
equal to 


= / u 


dx + v dy = 


of the compressible 



where the integration is extended over a closed curve around 
the profile P. By (15) 


Integrating over 


T c = a 0 Jn s - 
a closed curve 



&'(£)d£ 


around H. . 


But 


by (29) 





-4ttP sin cr 
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Introduce the values of R, O’, and 5 given "by (27), (28), 
and (40) to obtain the equation 

_ . */n s - 1 f . r -i \ n l net 

r„ = -4 tt" a. cos <! na + (n - 1)— > tan 

c 2-n° L 2 J 2-n 

or, by (16) 

F = -4Trq fS> — i cos -fna + (n - 1) — "l tan — — (47) 

u 2 ~ n 2 J 2 — n 

n is given in terms of M^, by formula (43). 

■ For an incompressible flow with the same velocity at 
infinity and the same position of a stagnation point the cir- 
culation is given by 

T ± = ~4rrq ej sin a, 

The values of r o / Fj, are given in table V and plotted in 

figure 7. It will be seen that compressibility results in a 
larger circulation. The additional circulation due to com- 
pressibility increases as the Mach number increases but it 
decreases as the angle of attack increases. 

The Kut ta-Joukowsky lift formula holds also for com- 
pressible flow. (See reference 7.) The lift (per unit span) 
is given by 

L " P«^ lr o f 

The dimensionless lift coefficient for a circle 


°L = 


1 2 
— p q r 

2 r co ^-co'. 


(r being the radius of the circle) is obtained from (47) 
and (16), noticing that this last formula holds for r = 1. 
Thus 


°L = 


8tt 


2-n 


cos 


{■ 


na + 


(n - 1) 


5 } 


tan 


na 


(48) 


2-n 



NACA TIT Ho. 970 


19 


Here a. is the angle of attack and a is determined by the 
stream Mach number. (See equation (43).) For an incompress- 
ible fluid 


Cjj = 8tt sin a. 


(49) 


The values of CL are given in table VI and plotted in fig- 
ure 8 . 


7. Comparison with the Karman-T sien Method 

The application of the linearised pressure-density rela- 
tion (Y ss -l) to compressible flow problems consists of 
two steps: (l) the solution of a boundary value problem for 

the equation of the minimal surface, and (2) the interpreta- 
tion of the results. The first step is purely mathematical. 
Concerning the second, the following point of view has been 
adopted in the present report. The velocity field (i.e., 
the values of q/a 0 ) computed under the assumption Y = -1 

Is considered as an apnr oximati on to the velocity field of a 
flow satisfying the actual equation of state (with Y = 1.405) 
which haB the same velocity at infinity. Accordingly, the 
stream Mach number H o is computed from q^/ag by means of 

the formula (9) obtained from the actual equation of state. 

Von Karman and Tsien (references 3 and 4) adopted a different 
point of view. They consider the fi'ctitious flow satisfying 
the linearized equation of state as an approximation to the 
flow of a compressible fluid possessing the same stream Mach 
number. Consequently, they compute the speed of the undis- 
turbed flow q # from the stream Mach number by means of the 
formula 



which follows from the linearized equation of state. (Fur- 
thermore, they take over from the results obtained by setting 
Y = -1 only the values of q/q .) Other interpretations of 
the results are also possible. Their relative merits can be 
determined only by comparison with rigorous solutions. If • 
equation (50) is compared with equation (9), it is seen that 
the theoretical and numerical results of the present report 
can be adapted to the Von Karman-T si en point of view by re- 
placing the value M^ wherever it occurs by the value 
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^co(Karman-T sien) 

CO 


(51) 


A * v + 1 

V S 


This relation is plotted in figure 9b. It is seen that for 
the values of the stream Mach number occurring in flows around 
a circular cylinder the difference between the two values of 
is quite small. 

Por the 
leads to the 

case of a circulation-free flow equation (22) 
following velocity correction formula 


q ^ <li 

1 - y, 3 


(52) 


oo ^-1 »eo ^ 



Here 

3 

y* = 

n - 1 
n + 1 

* 


so that , by (43 ) 





u = 



(53) 


/nzr^i + /TiZi^ 


q^ is the speed of the incompressible flow around the profile 

H (i.e. t of the flow with the complex potential G) . If 
profile distortion (the difference between the profiles n 
and- P) is neglected, q^ may be taken as the speed of an 

incompressible flow around P. Pormula (52) is the Von Karman- 
Tsien velocity correction formula, except that these authors 
obtain for y. the value 

■ ■ ■■ 

1 + /l - M ffl 3 

in accordance with their method of interpreting results ob- 
tained by setting V =■ -1. (The values at y. 2 are plotted 
in figure 9a.) 

Theoretically, formula (52) may be employed for circula- 
tion-free flows only. Nevertheless, It seems worth while to 
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try to use it for circulatory flows as well. Velocity dis- 
tributions computed by means of (52). are shown in figure 5. 
Figure 4 gives the values of the critical stream Mach number 
computed by means of (52). The agreement with results ob- 
tained by the method of this report is rather surprising. 

It is due to the fact that for circulatory flows around a 
circular profile, the critical values of the stream Mach num- 
ber are very low. Greater discrepancy should be expected in 
the case of slender profiles. 

Remark : The Glauert-Prandtl correction formula is de- 

rived under the assumption of a nearly uniform flow and 
should not be used for a circular profile. 


8. Comparison with the Method of Successive Approximations 


The method of successive approximations has been applied 
to the circulatory flow past a circular cylinder. The first 
approximation to the velocity potential has been computed by 
Lamb (reference 8), who used the Rayleigh- Janzen scheme and 
by Tamatiko and Umemato (reference 9), who used the Poggi 
method. In figure 5 velocity distributions computed accord- 
ing to the formula given by Tamatiko and Umemato are plotted. 
Figure 4 shows the values of the stream Mach number given by 
these authors. Tamatiko and Umemato also give a formula for 
the circulation 



1 + (ll. + i sin 8 a)M 2 
^12 3 / ® 


This formula leads to results quite different from those ob- 
tained in the present report. 


The second approximation to the circulatory flow around 
a circular cylinder has been computed recently by Heaslit. 
(See reference 10.) His numerical results are given in a 
form which does not permit an immediate comparison with the 
ones given in this report. 


GOHCLUDIHC- REMARKS 


It has been shown that the formula, transforming a cir- 
culatory incompressible flow around a closed profile into a 
compressible flow (obeying the linearized equation of state) 
may be used for the effective approximate computation of a 
purely subsonic flow around a profile closely approximating 
a given shape. 
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The following qualitative results have "been obtained 
for the case of a circular cylinder. It may he expected that 
they hold also for other shapes. 

1. The critical value of the stream Mach number M^ c 

decreases as the angle of attack a, (for the case of a cir- 
cle: the argument of the stagnation point) increases. 

However, J-3-M c /daf decreases as a increases. 

2. Compressibility results in a higher value of the cir- 
culation (and lift) than the one predicted by the theory of 
incompressible fluids for the same position of the stagnation 
points. This effect increases as the stream Mach number in- 
creases. However, it becomes less pronounced as the angle of 
attack increases, 

3. The Von Karman-Tsien velocity correction formula 
(which is theoretically applicable only for flows without 
circulation) yields good approximate results for circulatory 
flows of small stream Mach number and small angle of attack. 


Brown University, 

Providence, R. I., September 1, 1944* 
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APPENDIX 

THE INCOMPRESSIBLE FLOW AROUND THE PROFILE n 


The transformation 


z i 



+ e 


16 


- e 


-16 


(Al) 


takes the domain exterior to the profile n (see fig. 1) 
into a sector bounded by two rays through the origin which 
make the angles (tt/2 + 8) and -(tt/2 - 5) with the real 
axis. The transformation 


z a 


„ (tt/ 2 - 6) i 

Z x e 


(A2) 


rotates thi6 sector so that the lower ray coincides with the 
real axis. The transformation 
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e 21 VC s -a) 


C A3 ) 


maps the rotated sector into the' upper half-plane. The point 
£ = co is taken into Z L = 1 by the first transformation, 

into Z ss e^ 77 / 2 by the second and into Z = e^ by 

a 3 

the third, where 


X 


it/ 2 - S 
2 - n 


Finally the transformation 


Z 


He - lc r S 3 - 

Z 3 - .« 


(A4) 


(R, a real constants) takes the 
plane into the domain | Z| > R, 
the point £ = ») is taken into 
points, £ = e~^® and £ = -e^® 
and Z = -Re* 0 ”, respectively. 


upper half-plane of the Z 3 
The point Z 3 = e^ (i.e., 
Z = as . The two stagnation 
are taken into Z = Re“^ ff 


Next, the constants R and q- must be determined so 

that 


lim 

£ — > co 


lim /az 


dZg 

dZ x N 

<\d 0 £ -5> co\dZ 3 


dZ, 

4 0 


A simple computation yieldB the values 


(2 - ajcos a 2 

Thus, formulas (27) and (28). are verified. 

It follows that G( £) considered as a function of Z 
must have the form (3.5). 

In order to compute the correspondence .between the 
points of the profile H and those of the circle \Z\ = R, 
note that by (Al)t-o (A4) 
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j + e ls _ f Z + Re icr \ S ~ a 
l - e" 16 ^2 - He _lcr ^ 


(as) 


Equation (A5) is that of a Von Kar mkn-Tref ft a transformation, 
ae should have "been expected. Let £ = e^$ "be a point on 

the upper circle of the profile n and Z = Re* A the corre- 
sponding point on \Z\ = a. Substitute in (A5) and take ab- 
solute values: 


cos - 6 ) 

• 2 


cos i( A - <3") 
2 

sin i-($ + 6) 


sin — ( A + cr) 

2 


2 


(A6) 


This relation implies formula (32) for the angle i( A - c) . 

3 

Equation (A5) may *be written in the form 


+ Rj e 1 ^ fZ + Re icr ^ 

V.T _ -icrj 


s-n 


: i " H i 6 


-iP 


(A7) 


Z - Re 

where Rj_ and P are given by (25) and (26) and 

£i = £ + i(sin 6 - Ei sin p) 


Let £i = R x e 1 ^ be a point of the lower arc of the profile 

PI and Z = Re iA the corresponding point of the circle 
}Zf = R. Substituting in (A7) and taking absolute values 
yields 


cos i(\Jr - p) 

C 


cos — ( A - cr) 
2 

sin i(\J/ + P) 


sin i( A + cr) 

2 


2 


(A9) 


This relation implies formula (33) 


for the angle 


i(A- a). 


In order to compute the speed of the incompressible flow 
around fl, |G-‘(t)l. at a point of the profile, note that 
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For Z = Re iA , 


ls,(£)l \]1|M 

T 

- = 2 sin A + sir 

f 

4 

. a _ A. IT _ T>_i A 


(A10) 


n cr 


(All) 


Furthermore, for a point Z = Re corresponding to a point 
of the upper arc of fl (i.e., for -cr < A < rr + cr) 



dZ . R 

it 

dA 

t d$ 



dA 

c 0 s 6 

sin 

A + 

sin c 

d<$ 

(2 - n)cos cr 

sin 

® + 

sin 6 


and "by (A6) 


By (A10), (All), (A12), and (A13) 


(A12) 


(A13) 


|G«(e i$ )[ = 


2 cos 3 S (sin A + sin O') ' 
(2 - n) 2 cob 2 cr |sin $ + sin 81 


This is equation (30) . 

Similarly, for a point Z = Re* A corresponding to a 
point of the lower arc of n (i.e., for ct-tt<A< - cr) 

l4?l-#-l£l ui4> 


>f n 

(i. 

e . , 

for 

cr - 

TT 

dZ 

R 

dA 




it ” 

R i 

d\ir 




cos B 


sin 

A + 

sin 

cr 

n)cos 

~ 

sin 

V + 

sin 

P 


By ( A9 ) 


Thus, hy (A10), (All), (A14), and (A15) 
TJ « ^ ^ fi _ B ft \1 I _ j 1 


(A15) 


T.« . r» i\Jr . r . - _ . 6 \1 2 cos 3 B (sin A + sin <y) 3 

:G» R t e T - i ( sin 6 - R x sin B) 

J L J (2 - n) 3 cos 3 cr lain + sin B j 


This is formula (31). 



* 


fpATn-Jg II.- SOICTEI0NS OF BQPATIOS (36) 


n 

1.022 

1.024 

1.026 

1.028 

1.030 

1.032 

*.034 

1.036 

1.038 

1.040 

€ 

11.94° 

16.45° 

19.65° 

22.15° 

24.21° 

25*95° 

27.45° 

28. 7 8° 

29.95° 

31.01° 

r~ _ . .. 1 

n 

i.o4a 

1.044 

1.046 

1.048 

1.050 

1.052 

1.054 

I.O56 

1.058 

1.060 

€ 

31.96° 

32.82° 

33.62° 

34.34° 

35-01° 

35.63° 

36.21° 

36.75° 

37-25° 

37*72° 


n 

1.062 

1.064 

1.066 

1.068 

1.070 

1.072 

1.074 

I.O76 

I.078 

1.080 

€ 

38.16° 

38.57° 

38.96° 

39*33° 

39.67° 

4o.oo° 

40.31° 

4o.6o° 

40.88° 

41.15° 


to 
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TABLE III 


Radii of curvature of the profiles P 








*■ r 


TABLE I7.~ CRITICAL STREAM MACH MJHBER AS A FUNCTION 
OP THE POSITION OP THE STA®ATION POIHT 


a 

2° 

4° 

6° 

8° 

10° 

12° 

14° 

l6° 

18° 

20° 


0.402 

0.387 

0.373 

O.360 

0.348 

0.337 

0.327 

0.318 

O.309 

0.301 

- 

a 

22° 

24° • 

26° 

28° 


32 0 

34° 

36° 

38° 

4o° 

H .c 

0.293 

0.286 

0.279 

0.273 

0.267 

0.262 

0.257 

O.252 

0.248 

0.244 


a 

42° 

44° 

46° 

48° 

50° 

52° 

54° 

56° 

58° 

6o° 

H r 

OOi C 

0.240 

0.236 

0.232 

0.229 

0.226 

0.223 

0.221 

0.219 

0.217 

0.215 


w 

to 
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TABLE V 


Ratio of the circulation of the compressible flow 
to the circulation of the incompressible flow. 



.10 

.15 

,20 

.25 

.30 

.35 

1° 

1.01489 

1.03357 

1.05984 

1.09383 

1.13570 

1.18564 

20 

1.01475 

1.03325 

1.05926 

1.09290 

1.13431 

1.18366 

30 

1.01461 

1.03293 

1.05869 

1.09198 

1.13293 

1.18169 

40 

1.01448 

1.03262 

1.05812 

1.09106 

1.13156 

1.17975 

50 

1.01434 

1.03231 

1.05756 

1,09015 

1.13020 

1,17782 

74 ° 

1.01401 

1.03155 

1.05616 

1.08790 

1.12684 


100 

1.01367 

1.03079 

1.05477 

1.08567 

1.12352 


160 

1.01302 

1.02928 

1.05203 

1.08127 

1.11696 


SO® 

'1.01236 

1.02777 

1.04930 . 

1.07689 

1.11044 


30° 

1.01102 

1.02472 

1.04377 

1.06802 



40° 

1.00962 

1.02153 

1.03798 




45° 

1.00888 

1.01984 

1.03494 





TABLE VI 


Lift Coefficient 


\u 

.10 

,15 

.20 

.25 

.30 

.35 

<x\ 







i° 

.4452 

.4534 

.4649 

.4798 

.4982 

.5201 

2° 

.8901 

.9063 

.9291 

.9586 

.9949 

1.0582 

3° 

1.3346 

1.3587 

1.3925 

1.4363 

1.4902 

1.5543 

40 

1.7786 

1.8104 

1.8551 

1.9128 

1.9838 

2.0683 

5° 

2.2219 

£.2612 

2.3165 

2.3879 

2.4757 

2.5800 

7£° 

3,3264 

3.3840 

3.4647 

3.5688 

5.6966 


100 

4.4239 

4.4986 

4.6033 

4.7382 

4.9033 


150 

6.5895 

6.6953 

6.8433 

7.0335 

7.2656 


200 

8.7021 

8.8347 

9.0197 

9.2568 

9.5453 


30° 

12.7049 

12.8770 

15.1164 

13.4218 



40° 

16.3104 

16.5028 

16.7686 




45° 

17.9293 

18.1241 

18.3924 
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Figure 4.- Critical stream .Mach number as a funotion of the 
position of the stagnation point. 
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Figure ?a.- Circulation as a function of the position of the 
stagnation point. 
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Figure 8a..- Lift coefficient as a function of the position 
of the stagnation point. 
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Fig. 8b 






NACA TN No. 970 Fig. 9 







